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ABSTRACT 

The  original  design  of  the  plenum  and  bellmouth  contrac¬ 
tion  arrangement  of  a  subsonic  cascade  test  facility  did  not 
produce  sufficiently  uniform  flow  conditions  at  the  bellmouth 
exit  plane.  Pneumatic  measurements  revealed  sizeable  blade- 
to-blade  variations  in  velocity  and  flow  angle.  A  finite 
element  numerical  analysis  of  the  inviscid  flow  field  was 
carried  out  which  both  confirmed  the  need  for  and  guided  the 
design  of  a  modification  to  the  bellmouth  contraction  from 
the  plenum.  Following  the  modification,  in  which  the  original 
contraction  was  changed  to  two  two-dimensional  contractions 
in  series,  and  a  program  of  development  tests,  acceptably 
small  variations  in  velocity  and  flow  angle  were  measured 
at  the  bellmouth  exit  plane. 
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I.  INTRODUCTION 


In  1963  construction  was  completed  on  the  Naval  Post¬ 
graduate  School's  subsonic  cascade  test  facility.  Its  primary 
function  was  to  be  used  to  provide  information  applicable  to 
compressor  and  turbine  blade  design.  As  described  in  Ref.  1, 
several  unique  features  were  incorporated  in  the  design  of 
the  cascade  to  ensure  a  two-dimensional,  periodic  flow  at  the 
test  blading.  Previously  the  lack  of  good  two- dimens ionality 
in  other  facilities  had  plagued  efforts  to  correlate  air 
turning  angles  and  blade  surface  pressure  distribution  data 
[Ref.  2]. 

Initial  evaluations  of  the  facility  were  carried  out  and 
reported  in  references  1  and  3.  The  test  section  entrance 
flow  was  examined  and  initially  found  to  be  unacceptably  non- 
uniform.  Several  measures  were  taken  at  that  time  to  correct 
the  deficiency,  but  an  examination  of  subsequent  data  revealed 
that  only  partial  success  was  obtained.  In  particular,  an 
investigation  was  performed  using  NACA  65  series  blades  as 
inlet  guide  vanes  [Ref.  4] .  Difficulty  in  obtaining  two- 
dimensionality  was  reported  and  a  strong  recommendation  was 
made  that  the  test  section  inlet  flow  conditions  be  improved. 

In  preparing  for  tests  of  new  compressor  blade  designs, 
the  goal  of  the  present  study  was  to  take  whatever  steps 
were  necessary  to  first  obtain  a  two-dimensional  and  periodic 
flow  at  the  entrance  to  the  cascade  test  section. 
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For  purposes  of  this  report,  a  two-dimensional  and  periodic 
flow  will  be  considered  to  exist  when  the  dynamic  pressure 
variations  are  less  than  +_  II,  total  pressure  variations  are 
less  than  +_  .51,  and  flow  angles  are  within  _+  1.5  degrees 
of  their  centerline  value. 

Preliminary  pneumatic  measurements  indicated  that  the 
flow  from  the  plenum  itself  was  not  uniform.  A  numerical 
analysis  of  the  flow  field  through  the  plenum  chamber  and 
contraction  bellmouth  was  therefore  performed  using  the  finite 
element  method,  (FEM) .  Based  on  the  results  of  the  prelimin¬ 
ary  measurements  and  the  FEM  analysis,  plenum  redesign  and 
development  test  programs  were  carried  through.  The  final 
pneumatic  survey  results  showed  acceptable  levels  of  dynamic 
pressure,  total  pressure,  and  flow  angle  variations. 

This  report  documents  the  experimental  and  analytical 
steps  taken  to  obtain  a  suitable  plenum  geometry  for  the 
cascade  facility.  Following  a  description  of  the  cascade 
and  measurement  techniques  given  in  Section  II,  the  steps 
taken  in  the  development  of  a  successful  plenum  arrangement 
are  first  outlined  in  Section  III  and  then  described  in  the 
following  three  sections.  In  order  to  present  a  clear  account 
of  the  development,  details  of  each  section  of  the  work  have 
been  confined  to  appendices.  First,  the  FEM  program  is  given 
in  detail  in  Appendix  A.  Data  from  the  original  plenum  con¬ 
figuration  are  given  in  Appendix  B  and  the  plenum  development 
test  results  in  Appendix  C.  Finally,  detailed  drawings  of 
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the  construction  of  the  modification  to 
bellmouth  are  given  in  Appendix  D. 


the  plenum  and 
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II.  FACILITY  AND  MEASUREMENT  APPROACH 


A.  CASCADE  AND  INSTRUMENTATION 

A  description  of  the  Naval  Postgraduate  School's  Recti¬ 
linear  Cascade  as  it  was  originally  designed  is  given  in 
reference  1.  The  layout  of  the  complete  facility  is  shown 
in  Figure  1.  The  addition  of  two  fine  mesh  screens  at  the 
bellmouth  entrance  to  improve  flow  stability  is  the  only 
modification  to  the  plenum  which  is  reported  in  reference  1. 

A  follow-on  study  into  the  cascade  performance  was  conducted 
by  Bartocci  and  is  reported  in  reference  3.  Plenum  turning 
vanes  were  installed  by  Bartocci  to  direct  plenum  inlet  air 
towards  the  bellmouth  entrance  and  to  decrease  the  total 
pressure  fluctuations.  Multiple  views  of  the  plenum  chamber 
and  the  bellmouth  contraction  arrangement  as  it  was  at  the 
inception  of  the  present  work  are  shown  in  Figure  2. 

All  testing  in  the  present  study  was  conducted  with  the 
cascade  test  section  walls  vertical,  with  all  inlet  guide 
vanes  (Fig.  1)  and  test  blading  removed.  The  effects  of 
changes  in  plenum  geometry  were  evaluated  primarily  by  con¬ 
ducting  pneumatic  probe  surveys  at  the  exit  of  the  bellmouth, 
hereafter  referred  to  as.  the  test  plane.  Additional  informa¬ 
tion  was  obtained  using  tufts  for  surface  flow  visualization 
in  the  plenum  and  from  static  pressure  distributions  displayed 
on  water  columns.  A  Datametrics  Series  100  VT  Air  Velocity 
Meter  was  also  used  to  measure  velocities  in  the  plenum. 
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B.  PNEUMATIC  PROBE  SURVEYS 


1 .  Probe  and  Installation 

The  installation  for  probe  surveys  at  the  test  plane 
is  shown  in  Figure  3.  A  wooden  two-by-four  was  used  to 
traverse  the  probe  in  the  blade- to-blade  direction.  The 
coordinate  system  adopted  is  shown  in  Figure  4.  The  X-wise 
displacement  of  the  probe  was  manually  set  by  sliding  the 
bar.  The  displacement  in  the  Y  direction  and  rotation  of 
the  probe  were  adjusted  using  the  traversing  probe  mount  shown 
in  Figure  5,  which  was  attached  to  the  two-by-four  and  incor¬ 
porates  potentiometers  for  the  measurement  of  the  immersion 
(Y)  and  yaw  (g)  displacements. 

The  five  hole  United  Sensor  D- 1 2 5  probe  shown  in 
Figure  6  was  used  to  determine  total  and  static  pressure, 
and  yaw  and  pitch  angle.  A  U-tube  water  manometer  connected 
to  the  P 2  and  P^  ports  was  used  to  indicate  when  the  probe 
was  correctly  aligned  with  respect  to  yaw  angle.  An  approx¬ 
imate  value  for  the  static  pressure  was  obtained  by  pneumat¬ 
ically  averaging  ?2  and  P_.  An  exact  calibration  of  the 
probe  was  not  carried  out  since  only  relative  measurements 
(showing  the  degree  of  uniformity  across  the  test  plane)  were 
required . 

2 .  Data  Acquisition  and  Reduction 

Probe  and  other  reference  pressures  were  multiplexed 
via  a  Scanivalve  to  a  single  transducer  referenced  to  atmos¬ 
pheric  pressure.  The  transducer  and  Scanivalve  were  connected 
to  a  B  and  F  scanning  data  system.  Signal  conditioning 
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circuits  in  the  data  system  were  adjusted  so  that  pressures 
were  recorded  in  units  of  10  ^  inches  of  water.  Data  was 
printed  and  also  output  on  punched  paper  tape. 

A  Hewlett-Packard  Model  9830A  Programmable  Calculator 
System  with  a  two  Megabyte  mass  memory  disc  unit  and  X-Y 
plotter  were  used  to  process  the  probe  survey  data  from  the 
paper  tape. 

The  measured  quantities  and  uncertainty  intervals  are 
given  in  Table  I.  The  uncertainties  in  the  pneumatic  quanti¬ 
ties  were  calculated  from  the  2 a  intervals  of  25  consecutive 
readings  at  the  test  plane  center.  The  position  and  yaw 
angle  uncertainties  were  determined  by  the  accuracy  to  which 
the  appropriate  scales  could  be  read. 

Since  test  plane  Mach  number  did  not  exceed  .25, 
compressibility  effects  were  ignored.  Test  plane  Reynolds 
number  was  1.5  x  10^  per  foot  scale,  with  variations  of  +5%. 
In  addition  to  the  measured  yaw  angles,  the  following  (rela¬ 
tive)  dimensionless  flow  properties  were  evaluated  from  the 
measurements  in  Table  I,  and  their  distributions  across  the 
test  plane  were  examined: 

Dynamic  Pressure  (_Q) 

Q  =  ,  uncertainty  =  +_  .01, 

^REF 

where  Q(X,Y)  =  PT(X,Y)  -  PS(X,Y) 

Qref  =  PtC3Q’  5)  '  V30’  5)  =  Q(30’  5) 

and  the  uncertainty  in  Q  is  +  .008. 
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Total  Pressure  (AF^) 


_  APt(X,Y) 

APt  =  — j~ -  ,  uncertainty  =  +_  .05, 

1  ^REF 

where  APT(X,Y)  =  PpL  -  PT(X,Y)  , 

and  the  uncertainty  is  +_  .  049. 

Static  Pressure  (APg) 

_  apscx,y) 

APq  =  — ~ -  ,  uncertainty  -  +_  .015, 

5  ^REF 

where  APg(.X,Y)  =  Pg(X,Y)  -  P^  ,  uncertainty  =  +_  .01. 

Note  that  the  uncertainties  in  AP^,  and  APg  are  rather 
high  since  these  quantities  are  derived  from  small  differences 
between  two  larger  quantities.  Note  also  that  the  values  of 
QrEf  changed  over  the  period  of  a  test  run  due  to  variations 
in  the  plenum  pressure.  The  value  of  Qppp  at  each  (X,Y)  was 
therefore  derived  from  the  plenum  pressure  reading  taken  at 
the  time  of  the  Q(X,Y)  measurement.  This  technique  assumed 
a  linear  relationship  between  Q(30,  5)  and  Pp^> 

C.  OTHER  TECHNIQUES 

Qualitative  information  which  was  simple  to  obtain, 
greatly  accelerated  the  development  of  the  final  plenum  con¬ 
figuration. 

1.  Tufts 

Tufts  were  attached  to  the  plenum  walls  and  were 
observed  and  photographed  using  a  Polaroid  camera  during  testing. 
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(Average  velocities  of  less  than  10  ft/sec  allowed  the 
presence  of  an  observer  in  the  plenum.)  A  tufted  probe  was 
also  used  extensively. 

2 .  Static  Pressure  Distribution 

The  uniformity,  but  particularly  the  steadiness,  of 
the  cascade  sidewall  static  pressures  displayed  on  water 
manometer  columns  were  monitored  as  changes  in  configurations 
were  made. 

3 .  Velocity  Meter 

Velocities  from  the  sound  baffles  in  the  plenum  were 
measured  using  a  Datametrics  Series  100  VT  Air  Velocity  Meter. 
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III.  PROGRAM  OUTLINE 


The  goal  was  to  obtain  a  plenum  and  contraction  arrange¬ 
ment  which  would  give  at  the  test  plane  (outside  the  boundary 
layers)  maximum  peak-to-peak  variations  of  +  1%  in  dynamic 
pressure,  +_  . 5%  in  the  total  pressure,  and  +1.5  degrees  in 
the  flow  yaw  angle. 

The  program  was  carried  out  in  three  steps.  First,  tests 
were  made  to  determine  the  performance  of  the  existing  plenum 
configuration.  Second,  a  finite  element  analysis  was  carried 
out  which  confirmed  that  the  existing  bellmouth  could  not 
generate  a  sufficiently  uniform  flow  at  the  test  plane.  The 
analysis  also  guided  the  design  of  a  suitable  modification 
to  the  plenum.  Third,  an  experimental  program  was  carried 
out  during  which  a  major  modification  was  made  in  the  plenum 
geometry  and  a  plenum  and  contraction  configuration  which 
gave  satisfactory  flow  was  subsequently  developed. 

These  three  steps  in  the  program  are  documented  in  the 
following  three  sections  of  this  report. 

A  summary  of  the  configurations  tested  during  the  program 
is  given  in  Table  II.  Effects  of  configuration  changes  were 
evaluated  at  the  test  plane  by  conducting  blade-to-blade 
probe  surveys  ((first  at  50%  span)  and  spanwise  surveys  at 
selected  locations.  The  other  measurement  techniques  listed 
in  Section  III.C.  often  were  used  first  to  determine  whether 
a  particular  configuration  should  be  measured  in  detail. 
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IV.  PRELIMINARY  TEST  RESULTS 


Blade- to- blade  probe  surveys  were  obtained  with  the 
original  plenum  arrangement  (Configuration  0)  at  10,  30,  50, 
and  90%  span.  The  results  at  50%  span  are  shown  in  Figures 
7a,  7c,  and  7e.  (The  remaining  data  are  given  in  Appendix 
B.)  The  results  showed  unacceptable  variations  in  total  and 
dynamic  pressure  in  both  the  blade- to-blade  and  spanwise 
directions.  Spanwise  variation  of  the  various  parameters 
can  be  seen  by  comparing  values  at  given  blade- to-blade  loca¬ 
tions  in  the  results  at  10,  30,  50,  and  90%  spans.  The  yaw 
angles  were  also  outside  the  required  bounds. 

The  source  of  the  total  pressure  variation  was  traced  to 
the  screens  at  the  bellmouth  inlet.  Clearly,  the  normal  flow 
into  the  bellmouth  would  generate  a  velocity  which  would  vary 
in  magnitude  and  direction  over  the  face  of  the  screen, 
thereby  generating  a  corresponding  variation  in  the  total 
pressure  drop. 

The  screens  were  removed  (Configuration  1  )  and  the 
measurements  repeated.  Blade- to-blade  survey  results  are 
shown  in  Figures  7b,  7d,  and  7f.  A  spanwise  survey  was  con¬ 
ducted  at  X  =  30  inches  and  is  shown  in  Figure  8.  The  total 
pressure  variation  (_AP-p)  across  the  test  plane  was  seen  to 
have  decreased  to  +_  .5%.  The  variation  in  was  within 
tolerance  in  the  spanwise  direction;  however,  an  unaccept¬ 
able  gradient  still  existed  in  the  blade- to- blade  direction 
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at  all  spanwise  stations  surveyed.  The  yaw  angles  in  Figures 
7e  and  7f  suggested  tha.t  the  flow  was  still  converging  through 
the  center  501  of  the  test  plane.  In  addition,  observed 
spanwise  variations  in  flow  angle  near  the  ends  in  the  blade- 
to-blade  direction  suggested  the  presence  of  large  scale 
vorticity.  This  was  verified  using  a  small  string  tuft  at¬ 
tached  to  the  end  of  the  pneumatic  probe. 

The  question  then  was,  -  was  the  measured  velocity  field 
simply  the  result  of  the  geometry  of  the  bellmouth  contrac¬ 
tion,  or  could  something  be  done  in  the  plenum  itself  to 
correct  the  profile?  An  analysis  of  the  flow  from  the  plenum 
was  required  to  answer  this  question. 
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V.  FINITE  ELEMENT  FLOW  MODEL  PREDICTIONS 


A  two-dimensional  finite  element  analysis  was  used  to 
analyze  the  potential  flow  field  through  the  plenum  and  bell- 
mouth  contraction.  Two  questions  were  of  specific  interest. 
First,  what  was  the  velocity  profile  created  by  the  bellmouth 
contraction  at  its  exit  plane,  and  secondly,  what  distance 
downstream  was  required  to  attain  an  acceptable  uniform 
velocity  profile?  The  details  of  the  analysis  and  the  Fortran 
programs  developed  for  the  IBM  360  to  answer  these  questions 
are  given  in  Appendix  A. 

The  plenum,  bellmouth,  and  exit  passage  were  modeled  in 
two  ways  as  shown  in  Figure  9.  Model  1  was  the  most  simpli¬ 
fied  representation,  and  was  used  to  obtain  a  first  estimate 
of  the  velocity  profile  developed  by  the  bellmouth  contrac¬ 
tion.  Due  to  the  vertical  axis  of  symmetry  in  Model  1,  only 
half  of  the  model  needed  to  be  included  in  the  calculations. 

A  more  realistic  model  was  that  shown  in  Figure  9  as  Model 
2.  In  Model  2  the  flow  entered  the  plenum  from  one  side  as 
in  the  actual  plenum. 

The  computer  analysis  of  each  model  was  performed  in  a 
similar  manner.  Namely,  the  exit  boundary  condition,  which 
specified  flow  to  be  parallel  to  the  exit  walls,  was  applied 
at  boundary  planes  located  at  progressively  increasing  dis¬ 
tances  downstream  of  the  bellmouth  exit  plane  (Fig.  9).  For 
each  boundary  plane,  the  computed  velocity  profiles  across 
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the  bellmouth  exit  plane  and  the  boundary  plane  were  compared 
with  the  results  from  the  previous  run.  The  following  simi¬ 
larities  were  noted  in  the  results  of  the  two  models: 

1)  The  exit  boundary  plane  velocity  profiles  were  always 
more  nearly  uniform  than  the  bellmouth  exit  plane  velocity 
profiles . 

2)  As  the  exit  boundary  planes  were  moved  farther  down¬ 
stream,  the  exit  plane  velocity  profile  became  more  nearly 
uniform.  The  bellmouth  exit  plane  velocity  profiles  changed 
during  this  process. 

3)  As  the  exit  plane  velocity  profile  approached  a 
uniform  condition,  the  still  non-uniform  bellmouth  exit  plane 
velocity  profile  did  not  change  as  the  exit  boundary  was 
moved  even  farther  downstream. 

The  results  of  Model  1,  shown  in  Figure  10,  indicated 
that  at  exit  plane  (GH) ,  41.5  inches  downstream  of  bellmouth 
exit  plane  (CL),  the  velocity  profile  had  relaxed  to  give  a 
maximum  variation  of  0.3%.  The  velocity  profile  at  plane 
(CL)  had  a  variation  then  of  141  between  peak  velocity  at 
the  contour  wall  and  minimum  velocity  at  the  axis  of  symmetry. 
Similarly,  the  results  of  Model  2,  shown  in  Figure  11,  indi¬ 
cated  that  at  exit  plane  (LM) ,  50  inches  downstream  of  the 
bellmouth  exit  plane  (FS)  ,  there  existed  only  a  1,3%  velocity 
variation.  This  velocity  profile  was  not,  however,  symmetri¬ 
cal  due  to  the  asymmetry  of  the  flow  model.  Higher  velocities 
existed  on  the  inlet  side  of  the  exit  plane.  The  velocity 
profile  at  plane  (FS)  also  exhibited  the  asymmetry  and  a 
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15%  variation  between  maximum  and  minimum  velocities.  Figure 
12  shows  a  comparison  between  the  measured  profile  of  Q  for 
Configuration  1  at  50%  span,  and  that  calculated  from  the 
velocity  profile  of  Figure  11. 

While  the  magnitudes  are  different,  the  trends  in  the 
computed  and  measured  distributions  of  Q  are  clearly  similar. 
Since  the  analysis  is  a  two-dimensional  one,  the  magnitudes 
should  not  be  expected  to  fully  agree.  While  a  three-dimen¬ 
sional  inviscid  analysis  could  be  expected  to  be  more  repre¬ 
sentative,  the  inclusion  of  viscosity  would  be  necessary  to 
closely  approximate  the  actual  flow. 

There  is  sufficient  documentation  in  the  literature  to 
presume  significant  amounts  of  vorticity  exist  in  this  flow. 
The  rectangular  duct  shape  in  itself  can  create  sizeable 
corner  vorticity  [Ref.  5  and  Ref.  6].  Velocity  gradients  in 
the  boundary  layer,  under  the  influence  of  transverse  pressure 
gradients,  can  produce  sizeable  amounts  of  vorticity  also 
[Ref.  7].  As  was  stated  by  C.  Truesdell  [1954],  "Even  in 
nearly  irrotational  flows  the  relatively  small  amount  of  vor¬ 
ticity  present  may  be  of  central  importance  in  determining 
the  major  flow  characteristics." 

An  analysis  of  the  flow  in  the  orthogonal  plane  (i.e., 
in  the  spanwise  direction),  was  not  performed.  However,  an 
approximate  analytical  solution  for  a  similar  but  simpler 
model  sheds  light  on  the  probable  trends  in  the  velocity 
distribution.  For  this  analysis,  the  flow  field  is  created 
by  a  circular  cylinder  immersed  in  the  flow  between  two 
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parallel  walls.  This  is  illustrated  in  Figure  13.  The 
velocity  distribution  across  the  stream  at  the  plane  of 
maximum  contraction  is.  given  in  Ref.  8  as  the  following: 

sinh2(Trb/H)  cos  (2ttY/H) 
coshz  (ttX/H)  -  cos^  (ttY/H) 

1/2  s inh2  (irb/H)  sin“(2irY/H) 

[cosh2  (ttX/H)  -  cos2  (7tY/H)  ]  2 

The  important  result  to  note  is  that  as  the  ratio  b/H 
increases,  the  magnitude  of  the  velocity  variation  decreases. 
Using  a  similar  argument  for  the  cascade  bellmouth  contrac¬ 
tion,  one  would  expect  a  difference  in  the  maximum  velocity 
variations,  between  the  blade- to-blade  and  spanwise  direc¬ 
tions,  of  almost  an  order  of  magnitude.  This  was  consistent 
with  experimental  observation. 
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VI.  PLENUM  DESIGN  MODIFICATION 
AND  TEST  RESULTS 


The  measurements  reported  in  Section  IV  and  the  analysis 
reported  in  Section  V  suggested  that  the  bellmouth  contraction, 
as  designed,  could  produce  an  acceptably  uniform  velocity 
profile  in  the  spanwise  direction,  but  not  in  the  blade-to- 
blade  direction.  Several  perturbations  were  subsequently 
made  in  the  original  plenum  geometry,  and  the  test  results 
supported  this  conclusion.  A  complete  account  of  the  tests 
leading  to,  and  subsequently  verifying,  the  design  of  a 
modification  to  the  plenum  is  given  in  Appendix  C. 

The  basic  idea  in  the  modification  was  to  arrange  two 
two-dimensional  contractions  in  series,  each  with  sufficient 
length  to  allow  relaxation  to  uniform  velocity  of  the  accel¬ 
erated  flow.  The  spanwise  contraction  of  the  original  bell- 
mouth  contour  was  retained,  and  a  new  contraction  section  in 
the  blade- to- blade  plane  was  built  down  into  the  original 
plenum.  The  contours  and  construction  of  the  new  section  are 
given  in  Appendix  D.  The  final  configuration  (Configuration 
7)  of  the  plenum,  including  features  which  were  determined  to 
be  necessary  as  a  result  of  subsequent  development  testing 
(Appendix  C) ,  is  shown  in  Figure  14  and  Figure  15.  Referring 
to  these  figures,  the  baffles  were  installed  to  dampen  flow 
oscillations  observed  in  earlier  tests.  The  flow  diverter 
fastened  to  the  west  end  ceiling  was  installed  to  deflect 
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high  velocity  air  from  along  the  ceiling,  downward  through 
the  baffles  toward  the  plenum  center.  The  velocity  and  yaw 
angle  distributions  (Configuration  7  at  50%  span)  are  shown 
in  Figure  16.  Measurements  from  Configuration  1  are  shown 
for  comparison.  The  results  of  a  spanwise  traverse  at  X  =  50 
for  Configuration  7  are  shown  in  Figure  17.  The  variations 
in  both  dynamic  pressure  and  yaw  angle  were  seen  to  have  been 
brought  to  within  the  prescribed  limits  of  +_  1%  and  _+  1.5 
degrees  respectively. 
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VII.  CONCLUSIONS 


The  original  design  of  the  plenum  and  bellmouth  contrac¬ 
tion  arrangement  did  not  produce  sufficiently  uniform  flow 
conditions  at  the  bellmouth  exit  plane.  While  the  large  area 
contraction  ratio  C40:l)  ensured  uniformity  in  the  total 
pressure  when  no  screens  were  introduced,  the  original  bell¬ 
mouth  contraction  contours  in  the  blade- to-blade  plane  gave 
unacceptable  variations  in  velocity  and  flow  angle  at  the 
inlet  guide  vane  station.  The  inlet  guide  vanes  would  nor¬ 
mally  be  adjusted  to  control  air-inlet  angle  to  the  test 
cascade . 

Following  the  design  and  installation  of  a  modification 
to  the  plenum,  in  which  the  original  contraction  was  changed 
to  two  two-dimensional  contractions  in  series,  and  a  program 
of  development  tests,  acceptably  small  variations  in  velocity 
and  flow  angle  were  measured  at  the  inlet  guide  vane  station. 

Future  tests  will  involve  installation  of  the  inlet  guide 
vanes  and  an  evaluation  of  the  flow  downstream  at  the  cascade 
test  blading  location. 
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TABLE  I 

MEASURED  QUANTITIES 


Item 

Description 

Method 

Uncertainty 

X 

Blade- to-  Blade 
dimension 

X  =  0  in.  East  end 

X  =  60  in.  West  end 

Tape  measure 

+  .1  in. 

Y 

Spanwise  dimension 

Y  =  0  in.  South  wall 

Y  = 10  in.  North  wall 

Immersion 

Potentiometer  on 
probe  mount 

+  .01  in. 

3 

Flow  yaw  angle 

Yaw  angle 

Potentiometer  on 
probe  mount 

+  .25  deg. 

Pt(X,Y) 

Total  pressure  at 
the  test  plane 

Pressure  P^  on 

Pitot  probe 

+  .1  in.  H?0 

PPL 

Plenum  total  pressure 

Static  tap  in 
plenum  chamber 

V  ^  0 

+  .1  in.  H?0 

PS(X,Y) 

Static  pressure  at 
the  test  plane 

Pneumatic  average 
of  P?  and  P..  on 

Pitot  probe-5 

+  .  15  in.  Ho0 

pw 

Static  pressure  at 

X  =  30  in. ,  Y  =  0  in. 

Static  tap  on  south 
wall 

+  .1  in  1^0 

PATM 

Atmospheric  pressure 

Mercury  barometer 

+  .  01  in.  1^0 

PTARE 

Tare  reading  across 
pressure  transducer 

P  -  P 

a™  AIM 

+  .01  in.  H2O 

_ _ _ 1 

Note:  All  pressures  are  gauge  pressure. 
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TABLE  II 


PLENUM/BELLMOUTH  TEST  CONFIGURATIONS 


Configuration 

Number 

Description 

Figure 

Number 

0 

Original  design 

1  and  2 

1 

Removed  screens  from 
bellmouth  entrance 

1  and  2 
minus  screens 

2 

Simulated  "relaxation"  concept 
Added  plywood  vertical  wall 

C-l 

3 

Added  eastern  and  western 
lower  contractions 

C-4 

minus  fillets 

4 

Same  as  above,  but  moved 
turning  vanes  3  ft.  to  the  east 

As  above 

5 

Inserted  cardboard  fillets 
between  upper  and  lower 
contraction 

C-4  j 

? 

i 

6 

Removed  plenum  turning  vanes 

C-4 

minus  turning 
vanes 

7 

Installed  baffles  and 
flow  diverter 

14  and  15 

29 
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Figure  1.  Cascade  Facility. 
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Figure  2.  Plenum  -  Configuration  #0 
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Figure  3.  Bellmouth  and  test  plane 


in 

0 


33 


Figure  4.  Cascade  Nomenclature  (Test  Plane  at  Bellmouth  Exit) 
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Figure  5.  Pneumatic  Probe  and  Traversing  Probe  Mount 
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Figure  7a.  Blade  -  to- Blade  Dynamic  Pressure  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  230  ft/sec)  -  Configuration  #0 
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Figure  7d.  Blade- to-Blade  Total  Pressure  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  280  ft/sec)  -  Configuration  #1 
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Figure  7e.  Blade- to-Blade  Yaw  Angle  Distribution  at  the  Test  Plane 

(at. 50%  Span,  Velocity  =  230,  280  ft/sec)  -  Configuration  #0  and  #1 
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Figure  8b.  Spanwise  Total  Pressure  Distribution  at  the  Test  Plane 

(at  X  =  30  in. ,  Velocity  =  280  ft/sec)  -  Configuration  #1 
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Figure  8c.  Spanwise  Yaw  Angle  Distribution  at  the  Test  Plane 

(at  X  =  30  in. ,  Velocity  =  280  ft/sec)  -  Configuration  #1 
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Figure  9.  Plenum  and  Bellmouth  Finite  Element  Flow  Analysis  Models 
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Figure  10. 
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PNEUMATIC  DATA 
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Figure  12.  Comparison  of  Dynamic  Pressure  Distributions  (Pneumatic  Data 
at  501  Span,  Velocity  =  270  ft/sec)  -  Configuration  #1 
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Figure  14.  Plenum  -  Configuration  #7  (Side  Views) 


3  mh 
O  03 
co  pq 


50 


Figure  15.  Plenum  -  Configuration  #7  (Plan  View) 
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Figure  16b.  Blade- to-Blade  Yaw  Angle  Distribution  at  the  Test  Plane 

(at  50%  Span,  Velocity  =  270  ft/sec)  -  Configurations  #1  and  #7 
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Spamvise  Displacement  (incnes) 

Figure  17a.  Spanwise  Dynamic  Pressure  Disti'ibution  at  the  Test  Plane 
(at  X  =  50  in. ,  Velocity  =  270  ft/sec)  -  Configuration  #7 
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Figure  17h.  Spanwise  Yaw  Angle  Distribution  at  the  Test  Plane 

(at  X  =  50  in.,  Velocity  =  270  ft/sec)  -  Configuration  # 7 


APPENDIX  A 


FINITE  ELEMENT  FLOW  ANALYSIS 


A.l  Method 

A  computer  program  was  developed  to  solve  LaPlace' s 
Equation  (two-dimensional  potential  flow)  in  terms  of  the 
stream  function,  ip  .  The  purpose  of  the  analysis  was  to  sub¬ 
stantiate  the  premise  that  the  measured  velocity  distributions 
in  the  blade- to-blade  directions  were  a  result  of  the  inlet 
contraction  design.  The  finite  element  method  analysis,  FEM, 
was  used  to  solve  this  continuum  problem.  The  detailed  pro¬ 
cedures  for  an  FEM  analysis  are  presented  in  Ref.  9,  with 
numerous  specific  fluid  flow  applications  in  Ref.  8. 

An  FEM  analysis  results  in  a  piecewise  approximation  of 
the  field  variable  in  the  governing  flow  equation.  In  the 
present  case,  the  flow  is  governed  by  LaPlace' s  Equation, 

V  ip  =  0,  where  ip  is  the  field  variable.  The  solution  was 
sought  for  the  flow  through  a  two-dimensional  contraction 
with  the  same  physical  shape  as  the  bellmouth  of  the  NPS 
Rectilinear  Cascade.  Two  different  models  were  used  to 
represent  the  incoming  flow  and  plenum  boundaries. 

The  first  step  in  applying  the  FEM  was  to  divide  the 
continuum  into  a  discrete  number  of  elements.  Model  1  is 
shown  in  Fig.  A-l  with  the  appropriate  boundary  conditions 
specified.  Only  one-half  of  the  contraction  required  analy¬ 
sis  due  to  symmetry.  The  discretization  of  the  area  ABGH 
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was  carried  out  with  a  framework  of  graphically  estimated 
lines  of  constant  velocity  potential,  <j>  =  constant,  and 
streamlines  of  constant  *P .  This  method  was  both  accurate 
and  efficient  since  it  gave  smaller  elements  where  velocity 
gradients  were  high  (or  boundary  geometry  changes  occurred, 
as  in  the  contraction  contour),  and  larger  elements  where 
gradients  were  small. 

The  next  step  was  to  assign  nodes  to  the  element  vertices 
and  to  select  an  interpolation  function  to  represent  the  var¬ 
iation  in  4*  over  the  element  [Ref.  9,  Chap.  5].  A  linear 
two-dimensional  interpolation  polynomial  was  chosen. 

The  third  step  was  to  determine  the  matrix  equations 
expressing  the  properties  of  a  single  element.  Galerkin's 
method  of  weighted  residuals  was  chosen  to  accomplish  this 
task  [Ref.  9,  Chap.  4], 

The  fourth  step  required  an  assembly  of  the  sets  of  matrix 
equations  for  each  element  into  a  large  system  matrix.  The 
system  matrix  was  modified  at  this  point  to  introduce  the 
specified  boundary  conditions  [Ref.  9,  Chap.  3],  The  boundary 
conditions  were  specified  in  terms  of  node  values  of  the  field 
variable  of  its  first  derivative. 

Fifth,  the  combined  set  of  system  equations  were  solved 
simultaneously  for  the  nodal  values  of  the  field  variable. 

In  the  present  case,  a  computer  library  subroutine  called 
"LEQTIF"  was  used.  However,  any  equivalent  linear  equation 
solver  could  be  used  here. 
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Once  the  values  of  ip  at  each  of  the  nodes  were  deter¬ 
mined,  the  final  step  was  to  obtain  a  value  for  the  average 
velocity  between  any  two  adjacent  nodes.  Since  ip  was  assumed 
to  vary  linearly  between  nodes,  and  the  velocity  was  given 
by,  for  example, 


the  average  velocity  was  calculated  using 


U 


i>2  -  'P 
Y2  -  Y 


1 


1 


A. 2  Computer  Program 

A  listing  of  the  Fortran  IV  FEM  program  "LAPFEM,"  written 
for  the  NPS  IBM  360,  is  given  in  Table  A- I.  Table  A-II  is  a 
listing  of  terms  not  defined  in  the  program  comments.  The 
data  required  for  input  is  given  in  each  "READ"  statement. 

The  data  used  for  Model  1  with  (GH)  as  the  exit  plane  (see 
Fig.  9}  are  given  in  Table  A-III. 

Referring  to  Fig.  A-l,  and  using  for  example  the  element 
number  1,  the  correspondence  between  element  and  system 
coordinates  was  as  follows: 
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7 


1 


1 


2 


2 


3 


System  Element 


General  Element 


ICOR  (1,1)  =  7 
ICOR  (1,2)  =  1 
ICOR  (1,3)  =  2 

In  general,  the  rules  are  as  follows: 

1)  Correspond  all  elements  in  either  a  clockwise  or  an 
anti- clockwise  manner. 

2)  If  an  element  has  1  side  on  a  boundary,  then 
ICOR  (1,1)  and  ICOR  (1,2)  will  be  the  two  nodes 
of  that  side. 

3)  If  an  element  has  two  sides  on  boundaries,  as  does 
element  #1,  then  the  correspondance  proceeds  anti¬ 
clockwise  along  the  bounded  sides  first,  as  shown 
above . 

The  finite  element  construction  for  Model  2  is  shown  in 
Fig.  A- 2.  The  data  used  when  the  exit  plane  was  at  (LM)  are 
presented  in  Table  A-IV.  The  results  are  shown  in  Fig.  11. 
As  one  can  see,  a  larger  number  of  system  elements  are  re¬ 
quired  when  there  is  no  axis  of  symmetry. 

The  values  output  from  "LAPFEM"  are  as  follows: 

1)  Values  of  ip  at  each  node  point. 

2)  Velocities  between  specified  nodes. 


58 


OOOOOOOUOOOOOOOOOOOCOC'-’O^O^O^OCUCOOCOOOOOO^OOOOO 

—  (\i  co  st  in  no  °o  O'  o  <-*  cm  ro  ^  in  o  Is-  oo  O'  o  —  cm  m  *4*  in  h-  co  O'  o  —i  cm  ro  in  >o  r-  co  0s  o  ro  ^  in  vO  co 
o  o  c  o  o  o  o  o  o —4  —i  —  —  p-4  —  — •  r-s  *-4  cm  <m  cm  c\i  cm  eg  cm  cm  (\j  c\j  m  ro  co  ro  rn  ro  ro  m  ro  ro  m-  <r  >r  ^  «4-  ^  >r  >r 
ooooooooooooooooooooooc3ooo^cuooooyoooooooooooooo 
OOOOOOOOOOOOOOOOOOOOOOO 000^00000 0^00 OOOOOOOOOOOO 
(XQ-cxQ-CLO.q.Q-Q-O-Q-Q-Q-Q.CLO-Q-Q-Q-q.Q-q^Q-CLQ-QlQ-Q-^-O-D-Q-Q-Q-O-Q.Q-Q-ClU-q.Q.Q-ClQ-Q-q.CL 

<  <t  <r  <■  <<  <  <*  <c  **  <<<<<<<  <i  <t  <  <i  <i  <<*<<  <i  <i  <*  <  <t  <<<*  <<  <i  <:  <<i  <<  <r<  <  < 


LU 
X 


o  co 
Q  I 

UJ  *H I 


o 

CM 


a 

z 

— 

H 

CO 


2 

a 

o 

QP 

sx 


pq 

IX 

x 

< 


pq 

o 

PQ 

< 


K  X 

«w 

Z21” 

Vi 

— 

UJ 

LUO*- 

u 

UJ 

— 

— -J 

X  3 

O 

m 

muj 

UJK 

— 

a 

— 

(M2> 

— 1 2  co 

CD 

UJ 

in 

i— t  •* 

lu^H 

vU 

> 

z: 

m 

QU> 

OCQ.LU 

_J 

o 

>>C 

<  2: 

CL 

O^1 

_JQjUJ 

o 

aj 

O 

X  i _ 1 

— 

O 

-o 

• 

LUO  LU 

CM 

CM 

HH 

—  LLJ 

Z 

Zz"* 

1 

— 

l Hm 

LU 

<hmCJ 

1 

— 

CM—h 

LU 

t-*tO 

UJ 

in 

H  » 

3 

a:  •" 

UJ 

LU 

♦H 

— 

h 

1— <  V> 

LU 

C-J 

CM 

zo 

LU 

UJ 

LJ 

UJ 

<w 

>OJ 

—  «— 

OJ 

OC  *LJ 

UJ 

P- 

UJ  UJ 

• 

<t3o 

z 

— 

CD  J 

3  r5 

O  00 

UJ  UJ 

z 

o 

— 

-LU 

-j  —i 

•  UJ 

200  *a 

— 

in 

—  > 

<t  <  • 

CD— h 

2  — i_>C 

— * 

<C 

in  •* 

>  >o 

K 

OjZ^(2:< 

CM 

<M  — 

•  t 

>---• 

M  LUQlu 

OO 

w 

•— »o 

UO  *OCO 

XO 

J—O^TCC^tO 

m 

X 

—'—'CM 

•  >o> 

oo 

<ZUJ 

• 

— 

mow 

CD <  t<»— 

•x-j 

->*- '-J-J II 

• 

O 

— 

nJJ>UJ 

^OJwUJ(J<LU 

LJO0UJ— r 

o 

• 

incno-c 

Z  •  -j 

•OJ> 

LUX  — <IUJ 

:  • 

cu 

ZO*>  *J-CD 

Uj  ^  CO 

CM  Z  »*Z 

<1  *  J-OLJ 

LUt— 

«^a 

z 

>-^-H 

XCCO  •— i 

hiz 

UJ  CJ  -oc>- 

z: 

LLOin  •" 

X  3cdc^uj»-< 

OXtO  >x 

<i 

VO 

-eg— 

LU2^  »-H 

-J  3 

o 

X 

o 

— -^O 

ZZU>UIZ 

-JLU2-  >X 

X 

OC 

cotn^ru 

<  X 

o»  *+  * 

Q-Xlu  #a!<l 

LU 

LU 

— 

'Ol)v 

iriui 

— 

c  zu<o 

z 

rsj 

— 

C0C02Z 

-M-OJ  *Q 

—  1 

in 

^z>o 

•  UJ  »-H  <  1— 1 

hO< 

Z_ia3z 

’  # 

t 

Z 

IZ2U201UI 

oouJLL'  32 

3 

o 

(M> 

XI- 

LUX  hO 

z 

— 

O—  -CD 

OOUJLOUJVOUJO 

>UJQ^UJcCOO 

00 

in 

-o—»— 

LUXLUXUJX  00 

J  JCJ  J 

H* 

t- 

^Q-oinoc 

C^OhQ}-  ILUJ 

O^MJ-ZZ 

*  c>*zincMO 

LJ  LJ  U 

oo 

to  aozaj  lu 
Xyj  — <XS  21 
Xt-KC^SILUfOUJ 
<i~.Z— 33_J  |  —I 
C£2:>-<OujUJ(MLU 
lj-h— — ^ 
OX<LJ  OPQQC 

a:  ozo 

Cluj  a:  ll<u- 

xo 

IOHU-  CM  M- 


x>o 

^CQ  II 


*H  CO 
tl  It 


*  Wt0w^4t0 

* 

*to  -luo a 

*  >O0CCD 

*  oou.i<aj 

*  V^^MUJ 

*0  -3  -O 
*UL— -  — O 

*  cinzms: 

*■  XsOOCM 

*  SCOh-hH  • 

*  ZwOO^H 

■frOLUZU 

*  xcllucq 


o 

CD 

a 

Z 

— 

o 

CO 


o 

CD 

o 


O- 

z 
to 

^a^nzu.  o*z 
O  o  ot uz*  •* 
il  a.  n  D  *j 

OujOujP^LU_JLL‘^‘  UJ 

*33  X^O^Z 

*  *-j  •  _ I  •  _!.>  IT 

*  o<o<a<-  •  * 

*Z>Z>Z>UJO* 

*  iz* 

*  uuajOjQjUJLUH  * 
^ItItII  w* 


11  II 

<£ 

UJ 

VO  • 

o 

LUo 

£ 

a  • 

oco 

2 

z 

z 

X2 

d  z 

LUC 

o  Ul 

hZ 

I  S 

OOX 

>-UJ 

1— 

OOZ 

LU 

LL  U- 

CD 

oo 

• 

•  • 

CJO 

LU-> 

zz 

o 

—  » 

liiOh 

—  — 

UJC3Z2 

XX 

OO  C 

CMCJ  — 

OZ  *3 

win 

z  o 

Wm 

•  •—* 

o  •*  •*  •* 

•oco 

cd  in  in  • 

o  •  t- 

O-H-^ 

•  CDI— ^ 

^  -  it 

ID  UJX 

—  — 

ZJH 

o  • 

>ZX 

CD  II  HO 

X<OoO 

Z  • 

LU2^LU 

200  •  CD 

3r)aro 

•*H-  o 

<*  lu— *o 

OZ  th- 

ozoz 

CD  UJ  CD  XI 

OS  O 

UJLU  LULU 

ZLU>X 

xnx 

-0X0 

CL  LUO— * 

Z  Da 

LLU-U-U- 

uou-<i.>-h 

ooco 

ZOUO 


in 

Z  CM 

— **-<ooo 
— in  •*  •  •  • 

Oh^OOOOoO 
^sO  H  II  II  II  II  II  II 
— *-■ - 

inh 


_J  *XX 
LUOoO 

zz  * 

•  •  •* 
—  *LJO* 

-tx2z* 

-  * 
-  -* 


hhhh 


hhhh 
2ZZZ  — 

LULU  LULU  U  1 

QOQQ  — 
—  —  —  — *  ^ 

-*< 


LU Ul LULU*  — ' 
<3- 

•»•  ooooooooooco  *  <iq;^’>>>^ccd^z*-'q£  .  *-»  i^vjtoco*  ue. 


UJ 

UJ 

•JJ 

*  s>^.o 

*  -HHMH-<io*ajoooouo^ooixoinin»  o 

a 

X 

X 

*LJK  —  *— 

-a-  0—<*  QLU-LJCOCDCOH-I  — 

03  X  -‘—ir  -frX 

> 

>■ 

> 

*  o*— o  - 

^uzuu(qoUJ  * 

-  -If  00024 

»— 

H- 

►— 

*  -4  x 

*  cd>cd3>uj>5:Z* 

X  X-*  CDCDCD>* 

♦  *  — . 

— -1 

— 

4^  * 

*  ^auuuUu>* 

^O^O^-H- 

*  * 

*  * 

*  ^CCZcqZcDNZ^ 

4-  —  — i*-1—  * 

*  * 

# 

*  #  O 

r*-  —  *  cm 

*  * 

* 

r—  f-g 

OOOOOOOOOOO  OCOO(jUyUOUOu 

O^OOOo 

59 


OOOOOOOOOOOOWOOOOOOOOOOOOOCOOOOOOOUCOUOOOOOOOOOO 
o  ,-4  f>j  cn  4*  in  vO  oo  0s  o  — «m  rn  in  >o  n-  oo  O'  o  —  cnj  m  st  u '  ^  r*-  oo  O'  o  —  <m  m  in  >o  r*-  co  O'  o  —  cm  n >  >t  in  >c 

in  in  m  in  in  m  m  in  in  in  >o  >o  o  o  vO  >o  o  r-  r-  r-  r-  r-  r-  r-  c-  r-  oo  go  co  uj  oo  cc uj  co  uj  qo  o'  0s  0s  O'  O'  O'  0s 

0000000000000000000000000000000000<->0<->000  00000000 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO^OOOOCJOOOOOO 

<XCLQ_CLQ_0-Q_Q-Q-C^Q-<X^CLQ-0-a.Q-U_Q_CLQ_Q-Q_Q_CLQ_Q_CLQ_CLQ-Q_Q_CLCLQ-CLQ_Q-Q-Q-0-Q-Q-U_Q_CL 

<  <<  <  <r<r  <r  <t  <  <t  <r  ^  <t  <  <  <r  <r  <r  <r  <  <  <”<  <  **  <r  <  41  <*•  <  <r  <  <  <*  <  <r  <r<  <  <i  <r<r<x  <r  <  <  <i  <x 


CM 


CM 


0 

H 

—  — 

H 

LL 

•  m 

U. 

* 

<N 

UJ  •* 

•> 

« 

• 

0-^-4 

• 

O 

— i 

>- 11 

II 

m 

h  -J 

II 

a. 

— 

• 

•* 

H^UJ 

UJ 

* 

in 

0 

»  ♦— » 

UJ 

X 

X 

►— « 

00 

UJ 

►  *. 

O 

0 

00 

>r 

UJ 

►— • 

X  — 

-J 

< 

•» 

•H 

0 

w 

o_j 

<1 

> 

m 

•* 

0 

QC 

►UJ 

> 

» 

z 

» 

O 

•  ►— * 

> 

II 

z 

O 

CH'— 

> 

CC 

00 

z 

0 

»— 4 

UJ 

CC 

<r 

• 

UJ 

UJ 

r 

aa 

<r 

CJ 

> 

CJ 

UJ 

K 

— 

z> 

a 

: z : 

• 

0 

-* 

<r 

CM 

—  f— 

dC 

X 

UJ 

z 

I— 

<N 

0 

•■»■ « 

X 

0 

UJ 

• 

0 

O 

X  • 

u 

CO 

h- 

z 

cu 

y£ 

—1 

UJ 

O'  — 

CO 

UJ 

UJ 

IL 

►— » 

►m 

O 

UJ 

a 

> 

W 

•  » 

— z 

—  X 

2 

ULl 

m 

m 

cr: 

CMr-J 

ox 

o< 

00 

f- 

h- 

•» 

c 

II 

CDo 

cos 

CD  *~~i 

UJ 

Z 

X 

0 

Q- 

00 

(DOC 

CD 

< 

00 

in 

•— • 

Z  - 

Z< 

z^ 

*—4 

H 

*• 

»  — • 

►0 

-a 

a 

►  — 

— * 

•  <— 

rH» 

UJ 

UJ 

<1- 

•  m 

^-4 

x  •‘in 

II  ► 

II  * 

II  ^ 

h- 

CC 

02 

z  • 

•» 

y>jn 

—  |-»X 

— •  >"»x 

—  hX 

Z 

— <r 

UJ»  <  — 

O  — fM 

-J 

HJJ  •* 

0**  **in 

0—  •‘in 

O'—  ►in  — * 

— <2 

zoa 

•—•rn— • 

UJ 

-  —X 

X  o'"*  •* 

CDO-^  •* 

coo—  ►—  z 

Z3>W 

OQ-Z 

i-  *-o 

—‘—00 

Oco^m 

co^in 

ocoi-4inzj> 

> 

zuj 

M  00 

<i*~i  ► 

a:  •. 

Zo-H 

ZZ^M 

zuu 

•^1— • 

0  11  X 

CC 

Q-om 

•^O  * 

^zz  - 

►ZQ  •* 

+  CCr-lh- 

u_  <r  •* 

CO— CM 

0 

Z(jw 

0^  ►>•  0 

•— *  •*>  •  0  *—4  **>•  i 

<T  II  *-« 

O  Oh 

-J 

0 

0  M  ho  O  H  hO  O  II  H<J  r-4  11 

II 

hO 

O  II 

»  •* 

»— « 

•“-X 

O—  II  cc^o 

*—  II  CD  ||  O 

•—  If  cc  II  II  1-4 

*-»oo 

►O 

IUZm 

x— in 

UJ 

•» 

X  ►CO 

f-4  —  ►  OJ 

+  cn 

•— «  •» 

•*UJ 

— -J 

X 

»  W  t 

0 

—1 

'O—  ► 

►  • 

•  •>«-«  • 

•>  •  ►— 

—  I-.UJ 

J2.« 

•  >(VJ 

GU 

JJ 

►_jin 

n— hqo 

^-‘hOC 

I- 

w> 

<TUJh-» 

X  *^1 

<1 

LLJUU 

m  UJ  — 

•— •*— ‘***Z*-»'— » 

UJ 

in— 0 

h- 

mZ 

f“—  ► 

w'— W 

^wU 

^^0  '-'UJ 

<r 

0  00> 

r-J*— •  » 

•• 

Z  — X 

OUUffiUJOOZOUJOOQcDUJZZ 

>0>'-‘ 

0 

UJ  UJ  >•  •> 

•.wX 

UJ 

-J«-H 

UJ_JOO 

OCD>OOOCD>ZOOCD>00>> 

zjzx 

UJ 

5.  XOO^ 

•  xco 

O 

UJ  H 

ZUJ  ► 

—  OOhO 

ZOhQ 

OOmOZZ 

^  UJ  •  UJ 

X 

Of— 

Q-  ►r-J 

z 

hJ 

uin-inin 

m*— »-icd^2— 

mCC^ZmM 

—> 

O— M  U-~— 

z—  ► 

UJ 

hmUJ 

-J  — —  vO 

#0w^  .  0 

ww  «  O 

• 

► 

1  -  oxm 

•  '—in 

a 

uj^  ►cn 

^COhOo 
in-^moxin  •*  H  •  • 
rmu  ^  >J*cm 't  o  o 

*■♦-4  LPi  »*N  ♦‘I— »  II  II 

a:#  in '—'-'O'-*  *0*^  in— — 


HO^'tXO-'^inxO^'OX— *MOOXstN  -OCOXZ  0~x—  *  o 
LLUJrgfrif-HfMUJcvicn^cMLycvirn^cMr^^^fMcMp^^  NHHcgO^Hininox  a. 

O  —  •‘'V.  I  »«N  I  ►  Xf\j(M  —  **N  ►'—  —  U-  *‘SMOO(MHff)(NJ  00 

-*0  ►  ►  ►sO'-'O  *  ►'O'—  •*  ►w'O^hfNi^  ►h  UU 

cDin^^i— cD^in^H-LTsinh- w|— w<r  inf— so  w O *— 
<u^uj<z^^ulj<qv^uli<^<iu<uj2;  w<uj<rOinf---*f--'-*f-#  a:*  ^uJh^h*%*hOHH 
XZ00HSZQCHjZQQh2:QQ5:H2|-ti^CSH-20^<ULI<UJ<t4  0^QQ.<liJ<iu<lM^ 
Qt:^<r<N-*Q<^<i<i^Q:^<t<r»-^Q^<r<Acx»-4a.K-»cjin<ra:i-^cr:-JQ2:H-2:J—  x#  o*  02 

uu.uuujc^au.ujmQ^CLL  LuuJa:cJLUiiJLJ(XLjaiu-  h-uuqccc  <  qc  »— cr  ►-«  a:  aLUj;:>> 

U-H^cccc3u-HHa:Lr-3ix^,ai  cl  jLLua.ut-a.  -£ll-*  -t  qcll-^ul  ujouolju_cj> 

—  aiu. 

x 


rn 


CJ 


o 

CJ 

mo 

•—I  CNJ 


'O 


CM 


OQ 

CM 


O 

O 

ON 

rn*— * 


CO 


o 


■*  CX)-iOQCUJc£OQ0JO 
■*  ►LLISLL^LL 

*  *  X 

o  in  o*  nr  in  o  in 

cm  cnj  cn*.  -*  m  <r  >r  in 

*  *  o 

CM 

000 


60 


OOOOuOOCo^OOOOOOOOOOCJOOOOOOOOOOOOOOOO(-)OOuOuOOOOO 

od  ^  o  ^  ^  ^  so  o  <\j  co  ^ '  so  ^  oo  0s  o  on  ^  ^  p*- oo  C' ^  ^  -4*  sc  ^  oo  ^  o  i-i<\i  ro 

o>  O' O' o  ^  o  o  o  o  o  o  o  o  »-h  ^  »-4  ^  ^  ^  »-4  cn  ^  rsj  (\j  c\j  cv  eg  (nj  (\j  eg  n  >  co  co  m  ro  m  co  ro  co  ^  ^  >r  >r  >r 

O  O  »“i  *■“<  r“i  ^^r*H  *“"i  r*^  *“i  f“i  l*i  fH  r*H  r*H  *— ^  fH  i-H  f-H •“^  ^ 

oooooooooooooooooyoooooooouoooouoooootJoooooouooo 

Q.CLq_Q_0_Q_Q_Q_Q_O_q_Q_Q_Q_q_Q_q_U-CLq_Q_Q_Q_CLQ_q_Q_q_O-Q_q_Q_q^Q-Q-G_Q-O_Q_CLCLQ_Q_q_Q_Q_Q_Q_ 


<  <r<  <i  <  <  <<x  <*<x  <  <■  <r  <  <*■  <r  <  <<  <  <<  <r  <i  <r  **>  <  <r  <r  <  ** - 


oc 

UJ 


U 

< 

• 

UJ 

CO 

i—i 

oc 

•» 

< 

t— 

r*H 

UJ 

It 

3 

_J 

N-4 

X 

•* 

O 

LJ 

-*» 

♦» 

> 

M 

CL 

LL 

— * 

1 

a: 

2 

OO 

LL 

00  —* 

> 

2 

-J 

a 

2  CL 

OO 

00 

*m0 

-  2 

•» 

2 

> 

oc. 

^  00 

in 

•* 

o 

H  2 

i 

—4 

X 

LL 

LL 

-5 

U_ 

(Nl 

II 

— i* 

II 

•*  H 

II 

» 

—•♦—1 

o 

O 

3 

a 

—  II 

3 

a. 

•  ^ 

-J— 

o 

O 

<1 

UJ 

—*  “D  *-« 

< 

2 

2— in 

WM 

r- 

00 

*  •* 

00 

O  • 

UL^ 

U. 

< 

W-*— »  — * 

LL 

2 

Hw(\i 

OO  UJ 

o 

CD 

O 

-J 

'TQIL  w— _« 

O 

•> 

1 — LL  fH 

>-J 

K 

K 

DO 

On>  1  ^(VJ^ 

*-4 

-)  00  e? 

00  UJ 

2 

a. 

O  *^00  •  UL 

2 

» 

00 

_j>  •* 

1  > 

O 

□ 

O 

>-H 

a  co— 

O 

<1 

LU 

Ooox 

»• 

o 

o 

HH 

2  u 

Q-h-  QOO^o 

>•— » 

OO 

D5 

00  ►m 

o 

—  o 

— • 

1— 

<  CO 

2  *  w  t00>w0oo 

h~ 

> 

—1 

W  •> 

— CD 

>-*«-*  00 

>— * 

< 

*  a 

OOO  — ‘LL O  • 

Z> 

00 

< 

ujwin 

lls: 

OLJ 

w02 

21 

2 

20D^  O^CD  fH 

-J 

> 

X  M 

— ^00 

•2 

O  *2 

2 

CL 

*  X*  •* 

•*  «F>-  It  — * 

a 

u_ 

»■—.—  •> 

2  — 

O  - 

>o  •* 

> 

O 

LL 

^ oo  ii  it  in.  in 

in* 

*  4»00st 

UJ^-* 

OLUr-H-^ 

U 

LL 

< 

*<**  » 

-■»  tl  -5-^  O'Oos 

*  fr¬ 

*—4 

0  **m 

l  »-i>— i^-m 

•  tl  *-<CO  •  II  nh 

CD 

21 

*  »-*  #-»' 

•<  II  -0^*0  11 

• 

it*  Li? 

W 

*inx 

II  mmw^  II  *■ 

Ohw 

t|  Uhhw 

It 

• 

CL 

*c* 

W  CJI/)—  *•  t»  •‘^‘nO'-'nO 

M 

-fr  UJ*  — > *  •‘CvJnO^1 

Hh4ws«»(JJIL|^^ 

CO  u 

—  CD  O 

*•-1 

O 

>— • 

*  -J* 

*  -Owwl- 

|  ZZ^Zh* 

f  iJJ  ‘-*-1  ✓rfi-W  U-J  »  Tt 

2^Z<fZ'«* 


■* 


■  OCDO  ,</)w~'~w'"'lU<UJ*.  _*  *  UJ#  '->_UJ<2;<'vJ>>>>— l 

*  u-*  tnam*  u.«<u.i-3:(-#  *  -j*  x#oi«ii-z>ir\222zai(- 


w  m  u  ^ 

hZ4  *  _j 

*  * 

1-t  w  tl  ooooooooon^cc^*  — Ht* 

LLU  II  ^U.U  It  ->»  *  _J 

•fr  O*  LJ  II  UlL  h»>-»CIOq:*  * 

<** 

*i-h<X0CLJ  It  LJ  It  II  M  It  UJQC 

hGOU^ 

*  2^  Q 

CJ  LL  00  O )  0 )  00  0 ; 

3U-3*  H* 

0* 

it*  CL  a 

*  <1* 

*  * 

*• 

*  3U. 

^** 

*  * 

*  ♦ 

ir 

* 

w  *  *■ 

it*  <M<- 

>r 

^*  -w 

ir 

-»• 

0 

O 

***o 

*  •«■ 

O' CM 

*  it*  «* 

ir 

* 

at 

aoo 

*  * 

0 

O 

in 

*  *  -w* 

* 

* 

m 

'Tt^ 

^-***C0 

in 

in 

o> 

*  * 

in 

OUu  « 

□uuu 

Ooo 

UUO 

61 


OUOOOCOOOOOOOUCOOOOOOOOOOUOO^OOOOOOOOOoacOoOOCOO 

in  vo  r- co  a>  o  c\j  m  >4- u  ^  no  ^  o  f-<  ^  m  in  ^  ao  O' o  c\j  m  ^  in  h- co  o  c\j  rn  vo  oo  o  ^  eg 
>f  ^  >4-  m  in  in  in  m  u '  in  in  m  in  >o  so o  >o  so  so  "<j  \0  'U  *o  r*  h  r-  r-  h  r-  o-  o-  r-  e-  oo  co  oo  go  cu  qq  qo  go  oo  <x>  o'  cr  o' 


OOO^OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO^OOOOOO 

Q.Q_Q_Q_Q_Q_Q_Q_CLD~CLQ_Q-Q_Q_Q_Q_Q_Q-Q_Q_Q_Q-Q_Q_LLQ_Q_(i.Q_Q.Q-(XQ-Q-Q_Q-Q_Q_Ca-Q_Q_Q_Q-CLQ-Q_Q- 

<I<<<<<<<<1  <<  <  <<<<  <<<  <<<<  <  <<<<  <<1<<<1<<<1  <<<<<<<<<<<< 


co 


x 

CO 


• 

f 

fJ 

LJ 

— 

f 

— J 

CO 

• 

IL 

F 

eg 

#► 

•* 

in 

0 

m 

vU 

m 

r-* 

co 

CO 

F 

•» 

W 

X 

• 

<x 

CO 

LL 

• 

0 

fJ 

•J 

LJ 

F 

0 

F 

LU 

Z 

PI 

f— » 

• 

< 

*■* 

f 

< 

•* 

— 

O 

-J 

■» 

■» 

m 

a» 

UJ 

•> 

fH 

► 

CO 

LJ 

(NJ 

X 

UJ 

♦— 1 

Z 

— 0 

CO 

CO 

U 

U) 

•* 

wk 

> 

r- 

fH 

Z 

r— 

m 

f 

z 

F 

UJ 

Z 

m 

— 

LU 

f 

OJ 

UJ 

00 

•* 

in 

f-« 

CO 

Z 

*-“• 

UJ 

CO 

fH 

0 

CD 

O 

LJ 

Q 

f-« 

(NJ 

FF 

m 

CL 

O 

— 1 * 

•» 

LL 

F 

CO 

LL 

a: 

QL 

z 

— * 

■* 

X-* 

LL 

X 

LU 

LL 

< 

< 

* 

•-in 

UJ 

CO 

oc 

LU 

* 

* 

z 

00 

—0 

LJ 

— 1 

ex 

O 

— 

UJ 

>» 

UJ 

inro 

LJ 

f 

LJ 

LJ 

— 

— 

UJ 

1 

a 

fH  «« 

•» 

O 

> 

• 

• 

* 

LJ 

(\JZ 

LL 

eg 

cc 

LL 

C\J 

eg 

H 

“5 

z 

aj 

h- 

h- 

* 

* 

LU 

ILO 

O 

m 

z 

LU 

<_l 

* 

* 

CD 

>• 

z 

CO  ► 

Z 

*> 

LU 

z 

z 

— 

UJ 

< 

X 

z 

0 

< 

eg 

> 

N 

UJ 

coin 

CO 

LU 

UJ 

O 

O 

H 

«— 

:s 

LL 

-nO 

< 

_J 

O 

<w 

►— 

II 

—  COWL 

► 

UJ 

CO 

00 

O 

UJ 

D 

H 

* 

i 

h- 

l— 

03 

CD 

O 

—  CD 

< 

hUQ 

Z 

z 

Z 

z 

< 

< 

_J 

U-, 

N- 1 

CNJ>Z 

LU 

CD 

LU 

UJ 

UJ 

^oc 

CX 

^(X 

UJ 

00- 

-> 

LL 

z 

— 

»- 

z 

z 

+  < 

< 

*  < 

> 

O 

in  •. 

UJ 

•> 

CO 

UJ 

— 

UJ 

-** 

* 

—  * 

00  UJ 

*-«CX< 

—1 

X 

> 

C\J 

-J 

•— 

— 

1“ 

1  >LJ 

Z 

IMZ 

UJ 

NO 

CO 

LU 

CD 

LU 

cm  eg 

cn 

CSJfO 

M 

«*0 

O 

LL'-UOLU 

#»4— 

* 

*  0 

0 

*  LJ 

X 

X-5«J 

*— * 

«Z^J 

LU 

—  m 

JI 

LU 

LU 

*  * 

* 

** 

UJ 

fLU 

H 

Z 00  WHO. 

H 

Z  <  UJCO 

LJ 

K 

fH 

—  C\J 

• 

LLX> 

< 

o:>jjz 

a. 

< 

LU  UjZ  >— »— % 

<  — 

O 

ego 

•» 

00 

• 

z 

OoOlu<uo 

zo 

-J 

Ja: 

--j^c\JrO-jX 

♦  -J-ICD  + 

+ 

CD  * 

CL 

LL^-£.  »-«.*£ 

VI  • 

-J 

uj<i-huj  •. 

•*  'CLJeg^ujwrg 

Ol 

(Njmvjmx  o  o  •**—  *ozo  o 

v,  #win(\J  LL  *  LULL— f  II  _J 

'-o  ^uj  *  z<mz  11  o^—  < 

HfJ  It  •  *•  mZ>Omh  n  II  -3UJUJO 


-  II 

o  •*  •*. 

hxXUJIOmmm<  I 
►iv  ro  >-H  LU  ^  ^  w  o— ‘ 


I  II  I 


**  SoQ^^CviaDf^oOCD 
C0*^*~4C0* 

I  OLD  <  #<^<f-»<<c\J 


<LL— '-HZ 

►-'vm 

^“5  OD 

CX  CC  DC 

“5  X*— *X»-.-5cO 

w  » yj  w  UJ  w  w 

Z  «-_UJ<»-» 

*  <->* 

O-^  hZz 

**0 

OOO 

^Xi-»X»-h-5CD  || 

N  II  II  II  II  II  II  II 

QCf  w|-ZH 

*  *  oouj 

1#— • 

UJ<  •  O 

OWJo 

cm  •— 1  co  fh  eg  co 

OoO  J*-»CLZ 

*  *0:3:0. 

p-4llu'<h*h 

i^srcotn 

mhh 

^->>XXX  11  LU 

11  ^^HtMf-icoegeg 

u-HHiuacLJCia*  ^  cdz> 

^  co-£Z 

MOCO 

II  II  II 

If  11  II  If  II  II  OX 

C1C<<<<<<< 

j>3:lloz*  wjqi- 

0>”Lj>ljo 

cxcj-  a 

mOX 

l-Hf\|m^HCMroUJ< 

<t — 1 1 — *  — 1 — 1 p  — • 

F  | 

LU*  *000^ 

q^QVJUD 

SLL  fQ 

CD  CD  UJ  O  LJ  O  CD 

LUUJLLiUJLLiLLjJJ 

X 

*  *  *  fH 

-—J 

*  »  *  ^ 

CJO-fr 

O 

* 

* 

*  *  *  * 

wno* 

f-4 

* 

* 

* 

>r  co  eg 

***  * 

~<* 

* 

* 

* 

m  in  m 

*  * 

* 

* 

* 

* 

UCU  0 

00 

O 

O 

O 

62 


ooooooooooooooooooooooooooococo^oocoooooocaoooou 

m  m  no  r-  oo  cr  o  cm  co  ^  o  r-  co  a  o  *-4  cm  m  ^  m  >o  r-  go  &  o  ™  co  ,4-  in  so  r**  co  a  o  ^  cm  co  >r  in  'O  r- 00  a  o 
O'  a  ^  a  w'  v-n  o  o  o  o  o  o  o  o  o  o »-«  *-4  — ♦  ^  ^  *-i  — *  <m  «\i c vj  cm  cm  (\i  c\j cm  cm  cm  m  rn  co  r  >  rn  co  co  m  n  1  oo  m- 

^  —4  *-H  ^  —4  ^  (\|  (M  cvj  CM  ^  cvj  CM  C\J  C\l  OJ  CM  CM  CM  C\l  CM  (\4  CM  CM  Co  OJ  CM  CM  CM  CM  c\|  CM  f\j  CM  c\J  CM  CM  CM  CM  <N  <\i  (\l  CM  (\j <M  CM  CM 
OOOOOOOOOOOOOOOOOOOCM^OOOOOOOOOOQO^OOOOOOOOOOOOOO 
a  XQ,Q_aa.U-Q_Q„cXQ.Q_  0.0.0. 0-0.0- 0.0.  Q.Q.  00-000-0.0-0.  Q.  CX  CL  Q» Q- Q.  O.  tL  Q- CL  Q.  Q- Q.  a.  OO  q.  0 

r  <<<<1<<<<<<<<<<<C<<<<< 


00 

Z  00 

— • 

a  uj 

— 

*■» 

►M  O 

9  • 

• 

9 

h-  a 

CM  CM 

CVJ 

CVJ 

»-«  2 

**«■ 

O 

* 

* 

00 

O 

** 

* 

* 

3 

z  z 

h 

O 

a  uj 

— * 

HH 

— 

— • 

-I 

O  UJ 

Q 

— * 

— i 

-S 

2 

3 

o 

>  K 

3^ 

o 

3 

LU 

a:  uj 

o 

<  CD 

>> 

> 

> 

00 

o 

1  \ 

> 

1 

<x 

z  i- 

(X 

— 

3  Z 

►-H3 

< 

-5 

oo 

O  < 

i—3 

LJ 

N* 

H- » 

OD  1— 

w  w 

z 

OO 

>> 

3 

> 

> 

oo 

>  Z 

*w  ' 

O 

w 

2 

X  c 

oo  oo 

CD 

00 

oo 

o 

9 

O  LJ 

CD^J 

00 

CD 

CD 

»— • 

o 

3 

<r< 

cc 

LU 

< 

< 

t— 

9  9 

<  o 

< 

>o 

w 

w 

*— 1 

—rCD 

O  UJ 

w.'— 

* 

Io 

w 

o 

> 

s: 

4  4 

— 

02 

-h 

z 

oz 

CX  3 

A 

o 

— * 

LJ 

OJ2 

9  9 

• 

<z 

9 

9 

o 

LL2^J 

imcm 

CVJ 

UOJ 

CNJ 

CM 

UL2: 

o<i 

*  * 

* 

UJ 

* 

* 

V 

03 

COH-t 

*  * 

* 

ars 

in* 

* 

oc 

UJ 

QCt-UJ 

O 

o 

— 

ot- 

Osj-^ 

< 

2 

a*- CX 

in 

— **—  in 

m 

LL  LU 

in*-*— 

O 

o 

KQ< 

o 

CD 

— i  — 

>< 

oz 

3X 

CO 

0*5 

3 

o 

UJOZ 

o 

3iX  O 

00 

UM- 

0^3 

O 

CDQC 

>^> 

1- 

K 

CD 

02 

t— 

f— w 

CD 

a 

o 

— 1  — »  — 

x  x  — •— »  — 

< 

»— •< 

OX 

X 

II  x 

UJ>CC 

9  90  • 

it  •  *o  • 

W 

CCb- 

DC  1 

O  1 

UJ 

OCX 

CM  CM  O  CM 

— -^(MCMOCM 

4 

KOO 

o  — 

o— 

X 

H-O 

CX^O 

VN  v 

'v 

— 

<Z 

-o 

H 

X 

002: 

— ■  — i 

Hi  3  —  —*^“^ 

• 

XO 

^cn*- 

— % 

ao 

LL.  ^<4 

m-5H^ 

wwhh3H1^ 

CNJ 

O 

o 

^-»W 

QC 

3 

►^-5 

XX^3 

* 

K 

tLU  X  *2: 

•x  • 

Z2  O 

<i 

hOu 

w-w-O^ 

* 

ZQ 

OCw'0< 

OwvO 

<■<  UL 

0-1 

ZCD> 

OOUJO 

Oo0  2^UJZ 

UjUJ 

UJ  •OO'V* 

LLloO 

** 

UJ< 

UJ  O 

>>  •> 

CDCD>>  9> 

CO 

Z2: 

9— *CD<  12 

•  CD<Z 

•  .*z 

ZD 

ZZCD 

oo-*o 

<<00—0 

CD 

00 

—  _|<OCNJ< 

— <o<i 

CM  CM*  CC 

OO 

UJZ 

1  1  -J  1 

— ^  1  I  -J  1 

w 

—too 

^jlun-'Uj  +  'r 

_lwLU  + 

+  +*o 

PNIUJ 

*  «J<£LL 

vwUJv 

f^oo 

LUoO* 

ljj  •-«  w  2:  cm  eg  lu  w  2  m  n  cnj  rn  *  iu 

*UJZO 

*  *M-i* 

-w  *  * 

CNJOQ 

<•* 

>-i  w  ^  a  ^-4  ^  r-4  l-i  W  O  CNJ  CM  CNJ  cn  * 

4-00 

*  3  * 

cm  cm  ^rn  O  t—  H-  pu  cm 

«ooo 

>  * 

wajj-tNl<i<i<<w|_rg<x<<<*  -J 

H-4 

*  >UJO0* 

— ♦-4LUcM^»^^-H^UjnJ 

-Jw 

•  •  • 

U.OO* 

LUCU^*  J-J  JUUJQC* 

2. 

*  U_2LLJ* 

LUwi 

OOO* 

H- -■>—<  * 

Q->»C?C\JUJUJUJujQ.Or^)LUUJUJUJ*  cc 

o< 

3* 

UJUJ  >UJ 

OOOOUJLU>UJ 

II  It 

II  It  II  *  o  * 

II  it  II 

II  ></)CM  II 

II  II  II  *  UJ 

fNrvjrn*  z 

xo 

*  QO-J* 

II  II  t—  11  O  II  1 . .  II 

^m-tcgco*  oo 

km  ii  — i^cm •— *c\jh-  ii  — »rn 

•—a' 

*oz<* 

^  **  m  *—  cm  co  i-t  <n  ►_  m 

fTIrOU-LLU.*  SO* 

HHwCMU-i^*-<cM<\j»-<rouJCMmc\jrn*  lu 

IZ 

*  Z<^* 

LLU-wU. 

^CMLLlLwlL 

UJ"W 

wlU,H  II  <<l<l<wcm  it  <<<<*0 

ao 

*  * 

-JJU.-JC 

.J i i  — i in  — i 

.j  juuujajir 

LLf-H— J*--P t-J«J 

MU-J2U 

— '-J— 1* 

LJ 

*  UJ* 

QJUJ  LU O UJ  LU  UJ^  UJ 

LULU 

* 

a* 

LU<UJUJUJLUI— LLi^UjLUUJUJ*  UJ 

■fr  X* 

* 

N* 

*  X 

*  t~* 

'*r 

*• 

*  K- 

*  * 

* 

* 

* 

*  * 

O 

in 

* 

* 

* 

*  * 

nT 

*• 

—4 

—4 

uoouu 

UUUOuuUUUUuC 

63 


ooooooooooooocuooo 
^  rsj  ro  m  o  t-  °o  o  ^  <\j  m  m  o  r*-  oo 
>r  >t  >r  >r  ^  >r  ^  in  « '  in  in  m  m  in  in  in 

CSJ  C>J  C\i  rsJ  C\1  (NJ  <N  ^  CO  ^  CSJ  CM  C^J  Cvi  CN  (\l  CV 

ooooooooo^oooooooo 

aQ-Q- CL  a  CL  a  CL  CL  Q- CL  a  CL  CL  CL  CL  CL  (X 

<aL<m<<<<i<t<i  <:<<*<<<<<< 


U- 

o 

< 


t/) 

> 

l/> 

o 

K 

z 

N* 

UJ 

I— 

z 


»-4 

^(NJCOCM 

m 

m 

U- 

^p^^OsJ 

CM 

fO 

LL 

<C<<< 

<t 

< 

UJ 

I  —1  i  l 

—1 

_J 

o 

UJUJUJUJ 

UJ 

UJ 

o 

+  +  +  + 

4- 

+ 

^^•Hcvjcn 


h~i  H-  H-4  -3  ^  -5  h-4  -5  ic  UJ  LU  UJ 
+  +  + 

<C<<<<<<<^-3^ 
CO  C/>  t/>  to  to  GO  to  O0  ^ ^ 

- 

to  to  to  CO  CO  to  CO  to  to  U.  U-  LL 
|1  ||  II  M  II  II  ||  II  II  tOlOU) 


UJ 


UJ 


UJ 


UJ 


CO  *  *7  ^  1-4  ^  to  to  OO 

H4#  ^-0^-?  n  ii  M 

*  QC*  ******** 

*  to* 

*  »-•*  a: 


* 

*  t/)  to  to  to  to  to  to  to  to  to  to  4^:  h-o 

•pr 

ir 

■*  (/)V)i/)VJcov)totototototooa:uj 

IT 

* 

ir 

"W 

O 

7T 

*■  o 

O 

* 

■«■  in 

in 

r-4 

uUu 

TABLE  A- I I 


DEFINITION  OF  TERMS  NOT  DEFINED  IN 

THE  "LAPFEM"  PROGRAM  COMMENTS 

ITYPE 

Specifies  type  of  boundary  condition 

NEL 

Total  number  of  system  elements 

NSNP 

Total  number  of  system  nodal  points 

I IEL 

Number  of  element 

ICOR 

A  term  to  correspond  single  element  node 
numbers  (1,2,3),  with  system  node  numbers 
(1, 2, 5, 4, 5,  ...  ,  NEL) 

A 

System  matrix  of  coefficients  of  <J> 

U 

Column  vector  of  unknown  values  of  ^ 
^1^2’  *  •  •  ’^NSNp) 

F 

Column  vector  of  external  node  values 

SYSA(I,II)  System  A  matrix 

SYSF(I,II)  System  F  matrix 

ELA  I  II  Elemental  A  matrix 


ELF  II 

Elemental  F  matrix 

AREA 

Element  area 

I, II 

Matrix  row  and  column  respectively 

X ,  Y 

X  and  Y  locations  of  each  node 
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Figure  A-l.  FEM  Model  1 


o 

II 


> 

o 

o 

II  +-> 


> 


INI 


CXI 


35 


Figure  A- 2.  FEM  Model 


APPENDIX  B 


PRELIMINARY  TEST  PROBE  DATA 

The  following  figures  are  data  from  tests  of  Configura¬ 
tion  #0  and  Configuration  #1.  The  results  showed  unacceptable 
variations  in  total  and  dynamic  pressure  in  both  the  blade-to- 
blade  and  spanwise  directions.  The  yaw  angles  were  also 
outside  the  required  bounds. 
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Blade- to- Blade  Displacement  (inches) 
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Figure  B-3a.  Blade- to-Blade  Static  Pressure  Distribution  at  the  Test 
Plane  (at  101  Span,  Velocity  =  230  ft/sec)  - 
Configuration  #0 
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Figure  B-3b.  Blade- to- Blade  Static  Pressure  Distribution  at  the  Test  Plane 
(at  30%  Span,  Velocity  =  230  ft/sec)  -  Configuration  #0 
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Figure  B-3c.  Blade- to-Blade  Static  Pressure  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  230  ft/sec)  -  Configuration  #0 
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Figure  B-3d.  Blade  -  to- B1 ade  Static  Pressure  Distribution  at  the  Test  Plane 
(at  90%  Span,  Velocity  =  230  ft/sec)  -  Configuration  #0 
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Figure  B-4.  Blade- to- Blade  Yaw  Angle  Distribution  at  the  Test  Plane  (at 

10,  30,  50  and  901  Span,  Velocity  =  230  ft/sec)  -  Configuration  K0 
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Figure  B-5a.  Blade- to- Blade  Dynamic  Pressure  Distribution  at  the  Test  Plane 
(at  101  Span,  Velocity  =  270  ft/sec)  -  Configuration  #1 
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Figure  B-6.  Spanwise  Dynamic  Pressure  Distribution  at  the  Test  Plane 
(at  X  =  15  in.,  Velocity  =  280  ft/sec)  -  Configuration  #1 
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Figure  B-7.  Spanwise  Total  Pressure  Distribution  at  the  Test  Plane 

(at  X  =  15  and  30  in.,  Velocity  =  280  ft/sec)  -  Configuration  #1 
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Figure  B-8.  Blade  -  to- Blade  Yaw  Angle  Distribution  at  the  Test  Plane  (at  10,  30, 
50,  70  and  901  Span,  Velocity  =  270  ft/sec)  -  Configuration  #1 
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Figure  B-9.  Spanwise  Yaw  Angle  Distribution  at  the  Test  Plane 

(at  X  =  15  in.,  Velocity  =  280  ft/sec)  -  Configuration  # 1 


APPENDIX  C 


PLENUM  DEVELOPMENT  TEST  PROGRAM 

Several  alternative  designs  were  proposed  to  obtain  a 
uniform,  steady  flow  at  the  exit  plane  of  the  original  bell- 
mouth  contraction,  which  was  the  location  for  the  inlet  guide 
vanes  for  the  cascade  test  section.  Some  considerations 
governing  the  selection  of  the  approach  to  be  followed  were: 

1)  It  was  desirable  to  maintain  the  benefits  of  a  large 
settling  volume. 

2)  A  simple  construction  was  needed,  one  which  was  easy 
to  alter  if  necessary. 

3)  Minimum  boundary  layer  buildup  was  required  prior  to 
and  within  the  cascade  test  section. 

4)  Equal  boundary  layer  growth  on  opposite  sides  of  the 
proposed  design  was  required. 

A  modification  of  the  plenum  and  bellmouth  was  carried  out 
which  satisfied  all  of  the  above  requirements. 

As  shown  by  the  FEM  analysis  in  Appendix  A,  the  velocity 
distribution  created  by  an  East-West  contraction  required  an 
additional  length  of  non-converging  duct  in  order  to  relax 
(to  a  uniform  condition).  Configuration  2,  Fig.  C-l,  was 
tested  as  a  rough  simulation  of  the  "relaxation"  concept. 

The  arrangement  was  generated  by  using  the  existing  plenum 
turning  vanes  and  a  sheet  of  3/4  inch  plywood.  The  measured 
distribution  of  Q  shown  in  Fig.  C-2  showed  a  2%  variation 
through  most  of  the  tunnel  width  except  in  the  deadwater 
regions  (evidenced  by  P^.  readings)  at  either  end.  The  yaw 
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angles  in  Fig.  C-3  were  similar  to  Configuration  1  except  as 
expected  at  the  ends  of  the  test  plane.  A  major  modification 
was  made  to  the  plenum  to  obtain  two  separate  contractions 
in  series  with  adequate  lengths  for  relaxation.  In  Configur¬ 
ation  3,  the  east  and  west  contours  of  the  bellmouth  were 
covered  by  vertical  walls  dropping  down  to  a  newly  constructed 
set  of  East-West  contraction  contours.  The  construction 
sketched  in  Fig.  C-4  (and  detailed  in  Appendix  D)  allowed 
some  relaxation  of  the  velocity  distribution  created  by  the 
new  East-West  contours  before  the  flow  entered  the  original 
North-South  bellmouth  contours.  The  cardboard  fillets  shown 
in  Fig.  C-4  were  not  installed  in  Configuration  3.  Measure¬ 
ments  of  Q  and  yaw  angles  were  made  and  the  results  are  shown 
in  Fig.  C-5.  The  non-uniformities  were  acceptable  through 
the  center  40  inches  of  the  test  section.  The  velocity  dis¬ 
tributions,  however,  peaked  rapidly  at  both  ends  and  the 
sudden  increase  in  yaw  angles  indicated  a  strong  influx  of 
fluid  from  the  East  and  West  walls  toward  the  tunnel  center. 
The  close  confinement  of  the  deviant  yaw  angles  to  the  tunnel 
ends,  at  X  =  10  inches  and  X  =  50  inches,  was  thought  to  be 
partially  caused  by  the  proximity  of  the  plenum  turning  vanes 
to  the  lower  contraction. 

In  Configuration  4,  the  vanes  were  moved  three  feet 
towards  the  east  end  of  the  plenum  and  the  measurements  were 
repeated.  The  results  shown  in  Figures  C-7  and  C-8  were 
similar  to  those  from  Configuration  3;  however,  the  peak 
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velocities  at  the  ends  were  decreased  slightly  and  the  yaw 
angles  varied  across  the  entire  40  inch  channel  width. 

At  this  point  the  west  end  lower  contraction  surface  and 
vertical  wall  was  tufted  with  three-inch  strips  of  light 
string  in  order  to  examine  the  surface  flow.  A  tufted  rod 
was  used  to  probe  the  flow  outside  the  wall  boundary  layer. 
Figure  C-9  shows  a  sketch  of  the  tuft  positions  with  air  flow 
through  the  test  plane  at  a  dynamic  pressure  of  15  inches  of 
water.  The  positions  were  traced  from  a  Polaroid  photograph 
taken  from  within  the  lighted  plenum.  The  view  was  taken 
from  the  east  end  of  the  plenum  looking  west  at  a  45  degree 
upwards  angle.  The  following  observations  were  made. 

1)  Flow  on  the  vertical  wall  in  the  vicinity  of  the 
upper  contraction  converged  more  rapidly  than  the 
contraction  walls. 

2)  Large  regions  of  deadwater  existed  in  the  upper 
corners . 

3)  Strong  corner  vortices  (verified  by  the  hand  held 
probe)  were  present  at  the  junction  of  the  ceiling, 
upper  contraction  walls,  and  vertical  walls. 

These  observations  pointed  to  the  presence  of  substantial 
secondary  flow  of  the  boundary  layer  fluid.  Cardboard 
fillets  were  installed  in  Fig.  C-4  to  eliminate  the  dead- 
water  areas.  Pneumatic  measurements  taken  with  this  con¬ 
figuration  are  shown  in  Figures  C-10  and  C-ll.  The  results 
indicated  little  change  in  the  Q  distribution;  however,  the 
peak  yaw  angles  were  notably  increased. 

The  effect  of  the  position  of  the  turning  vanes  was  ques¬ 
tioned.  When  the  plenum  turning  vanes  were  displaced  three 
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feet  to  the  east,  it  was  apparent  that  very  little  flow  was 
passing  through  any  hut  the  top  turning  vane  passage.  The 
vanes  were  therefore  removed  altogether  as  their  functional 
benefit  was  questionable.  The  results  of  probe  measurements 
for  this  Configuration  6  are  shown  in  Figures  C-12  and  C-13. 
The  distribution  of  Q  was  seen  to  be  similar  to  the  previous 
one  but  an  increase  was  noted  in  the  peak  values  of  the  yaw 
angles.  The  yaw  angles  measured  at  20%  and  70%  span  were 
consistant  with  the  proposed  existance  of  strong  contra¬ 
rotating  vortices  centered  at  about  X  =  10  and  45  inches. 

An  audible  low  frequency  pressure  pulsation  was  noted 
and  the  static  pressures  in  the  cascade  test  section  gave 
an  unsteady  measurement  of  +_  0.75  inches  of  water  on  the 
manometer  columns.  These  observations  suggested  that  the 
plenum  turning  vanes  had  served  as  a  damping  structure  and 
perhaps  as  a  spoiler  which  deterred  the  formation  of  large 
scale  vortex  motions. 

The  Datametrics  Air  Velocity  Meter  was  used  to  probe  the 
plenum  chamber.  Flow  velocities  through  the  sound  baffles 
were  quite  uniform.  However,  at  a  distance  of  eight  feet 
downstream  of  the  sound  baffles  approximately  75%  of  the 
flow  was  seen  to  have  passed  through  the  upper  one-half  of 
the  plenum  at  this  location.  It  was  also  noted  that  little 
flow  was  actually  moving  around  the  east  end  lower  contrac¬ 
tion  contour. 

It  was  clear  at  this  point  that  in  order  for  the  lower 
contraction  to  work  properly,  it  was  necessary  to  prevent 
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the  flow  from  forming  a  smoothly  accelerating  jet  from  the 
distributed  source  represented  by  the  sound  baffles,  to  the 
line  sink  represented  by  the  bellmouth  exit  plane. 

Consequently,  the  plenum  was  reconfigured  as  shown  in 
Section  VI,  Figures  14  and  15.  Baffles  were  positioned  to 
behave  as  spoilers  and  to  better  distribute  the  flow  below 
the  contraction.  A  4-by-12  foot  construction  of  1/8  inch 
masonite  was  installed  on  the  west  end  ceiling  to  divert  the 
high  velocity  fluid  from  the  ceiling  down  through  the  baffles 
The  measurements  of  Q  and  flow  angles  for  this  final  Configur 
ation  7  are  presented  in  Section  VI. 
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Figure  C-l.  Plenum  -  Configuration  #2 
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Figure  C-3.  Blade- to- Blade  Yaw  Angle  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  260  ft/sec)  -  Configuration  #2. 
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Figure  06.  Blade- to- Blade  Yaw  Angle  Distribution  at  the  Test  Plane 
Tat  50%  Span,  Velocity  =  250  ft/sec)  -  Configuration  #3 
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Figure  C-8.  Blade- to- Blade  Yaw  Angle  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  250  ft/sec)  -  Configuration  #4 
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Figure  C-9.  Tuft  Positions  on  West  End  Lower  Contraction 

(Viewing  East  to  West,  45  Degrees  Upward  Angle) 
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Blade- to- Blade  Displacement  (inches) 

Figure  C-10.  Blade- to- Blade  Dynamic  Pressure  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  270  ft/sec)  -  Configuration  #5 
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Figure  C-ll.  Blade- to- Blade  Yaw  Angle  Distribution  at  the  Test  Plane 
(at  50%  Span,  Velocity  =  270  ft/sec)  -  Configuration  #5 
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Figure  013.  Blade- to-Blade  Yaw  Angle  Distribution  at  the  Test  Plane 
(at  20,  50  and  70^  Span,  Velocity  =  250  ft/sec)  - 
Configuration  if 6 


APPENDIX  D 


PLENUM  MODIFICATION  DRAWINGS 

The  following  drawings  show  the  construction  of  the  east 
and  west  lower  contractions. 

Construction  Materials: 

Figures  D-l  through  D-4 

1.  All  framing  material  is  construction  grade  pine. 

2.  The  surface  covering  is  1/8  inch  untempered  masonite 
paneling  in  4  x  8  ft  sheets. 

3.  Countersunk  wood  screws  secure  the  masonite  to  the 
framing . 

4.  Framing  is  secured  to  the  cement  wall  with  power- 
driven  nails. 

Figure  D- 5 

The  upper  contour  is  part  of  the  original  bellmouth 
constructed  from  welded  1/8  inch  steel  plate. 
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Figure  D-2.  W-E  Elevation,  Western  Lower  Contraction  Framework 

Side  View  from  South  Scale  1:10 
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Figure  D-3.  E-IV  Elevation,  Eastern  Lower  Contraction  Framework 

Side  View  from  North  Scale  1:10 
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Figure  D-4.  Lower  Contraction  Contour 
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